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We report on a study of the evolution of magnetic fields of neutron stars, driven 
by the expulsion of magnetic flux out of the proton superconducting core of the 
star. The rate of expulsion, or equivalently the velocity of outward motion of flux- 
carrying proton-vortices is determined from a solution of their equation of motion. 
A determination of the effective forces on the fluxoids moving through the quantum 
liquid interior of neutron stars is however confronted with many ambiguities about 
the properties of this special case of superconductivity in the nature. Also, the 
behaviour of the fluxoids at the core boundary, and the subsequent evolution of 
the expelled flux within the highly conductive surrounding crust, are other related 
issues that have not been so far explored in any great details. 

1 Introduction 

Observational evidence has accumulated_that magnetic fields of neutron stars decay over time 
scales of tens of million years or more lJ'cB. Theoretical modelling of such a field decay calls 
for a study of the expulsion of magnetic flux out of the superconducting core of the star, which 
is embedded within the highly conductive crust of the star. The study would bear on an un- 
derstanding of the magnetic and superconductivity properties of matter under conditions much 
different than is usually realised for the terrestrial superconductors. Some of the novelties might 
be Hsted: 

• type I/II proton superconductivity; normal electrons, 

• "field cooling" in an intrinsic field, originally supported by electron currents, 

• coherent scattering of electrons off fluxoids, 

• "transfer" of flux from proton super-currents to electron currents outside, 

• freezing of fluxoids end points; geometry of the expelled flux, and 

• spherical geometry and km size of the superconductor and fluxoid lattice. 

The velocity of outward motion of flux-carrying proton-vortices (fluxpids.) may be determined 
from a solution of the Magnus equation of motion for these vortices llSES, taking into account 
the various forces that act on them. For the terrestrial type-II supercunductors flux movement 
is usually driven by the Lorentz force due to an applied transport current. In contrast, in the 
interior of a neutron star there are other forces acting on the fluxoids. These include i) a force 
due to their pinning interaction with the moving neutron vortices, ii) viscous drag force due to 
magnetic scattering of electrons. Hi) buoyancy force, and iv) curvature force. The derived radial 
velocity of the fluxoids at the core-crust boundary would in turn determine the rate of the flux 
expulsion out of the core. The magnetic evolution of the star may be thus determined by further 
following the Ohmic dissipation of the expelled flux within the crust. A detailed description 
of the modelling of such a flux expulsion and field evolutiniL pad the astrophysical implications 
of the results of our computations is reported elsewhere E3^Ei]^E3. Here we give a brief outline of 
aspects of the models directly relevant to the superconductivity effects. Our aim would be to also 
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highlight the basic issues, listed above, which require further theoretical investigations pertaining 
to the superconductivity and magnetism in neutron stars. 

2 The Superconductor 

All theoretical models of neutron stars predict a (quantum) liquid interior ("core") at densities 
(^ 2.4 X 10^"* g cm^'^) similar and above that of the nuclear matter, consisting of bare neutrons 
plus an admixture of protons and electrons, which is ^uixminded by a solid metallic crust of 
neutron-rich nuclei and relativistic degenerate electronsEl'Ej'EiJ. Neutron stars are "cold" objects, 
having internal temperatures T ~ 10® (K) << TFcrmi ~ 10^^ K. For the same theoretical 
reasons that terrestrial matter is argued to become superconducting with a transition temperature 
T-jiij 10^'^ TFormi the ucutrou stars are predicted to have superfluid-superconductor interiors 
Bla^ESl. Cooper pairs, in a neutron star, are however formed due to the long-range attractive part 
of nucleon-nucleon interaction among protons, as well as neutrons. The proton component of 
the core is usually assumed to form a type-II superconductor, based on its estimated values of 
the penetration depth Ap ^ 10^ fm and the coherence length ^ few fm. The possibility of 
type-I proton superconductivity, ^p > %/2Ap, has been also considered, allowing for the existing 
uncertaip|i|ti£S in the values of density, effective mass of protons, and transition temperature in 
the coreB'EZl. However, no attempts have been reported on the modelling of the field evolution 
of the star for the case of type-I superconductivity. The transition to the superconductirig state 
is believed to occur at a constant field, originally present in the core in its normal stateu. This 
should lead to formation of a lattice of proton flux tubes (fluxoids) which extend some few km 
across the spherical core of the star. Distorted geometries of lines might be also realized due 
to the presence of field components other than a pure dipolar field in the normal stateu, which 
is neglected in our models for simplicity. An interesting feature of this flux trapping is that, in 
contrast to the usual requirement for the Meissner effect, the field need not be expelled in spite 
of being smaller than the associated lower critical value Hd ~ 10^^ G, which is indeed the case 
for the typical field strengths of neutron stars ^ 10^^ G. The large conductivity of the matter 
in the normal state implies dissipation time scales, for the electron currents in the core, many 
orders of magnitudes larger than the time taken for the star to cool down below its transition 
temperature Tc. This is arEued to imply freezing in of the flux even at field strengths smaller 
than the lower critical value Efliil. 

In the thin layer (~ 1 km) of the crust, surrounding the superconductor core, matter at the 
density range 7 x 10^ (g cm"'^) < p < 2.4 x 10^"' (g cm^'') is highly conducting. The estimated 
Ohmic dissipation time for the crust, due to electron scattering by phonons and lattice impurities 
is uncertain and lies in the range 10^-10^ yr. Moreover, the unknown geometry of the expelled 
flux and itSjiransport behavior within the crust further complicates the calculation of its decay 
time scaleli3l3. The distinctive features of the flux expulsion in the present case, in contrast to that 
of the terrestrial superconductors, should be born in mind. Firstly, this is not an expulsion of the 
flux associated with an externally applied field that might be understood in terms of interaction 
of the fl-uxoid current with the London supercurrent at the surface, and a readjustment of the 
latter. The flux expulsion out of the core of neutron stars requires electron currents to be induced 
in the surrounding crust, at the cost of proton supercurrents of the fluxoids being annihilated. 
Furthermore, the fluxoids end points spout out in a highly conductive surrounding medium. 
Considering the spherical geometry of the superconductor, it is not clear whether the expulsion 
of the flux is realised only at the magnetic equator or elsewhere as well. Also, the back reaction 
of the flux accumulated in the crust, as well as the effect of a surface "barrier", on the fluxoids 
motion close to the boundary might have significant new results. For example, the surface barrier 
for the present case, having external fields smaller than Hd, shpuld be an outward declining hill 
causing the fluxoids to be drained out of the superconductor 13. Such an expulsion mechanism 
due to surface effects has not been so far addressed, in the case of neutron stars; one should see 
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how the corresponding rate of flux expulsion, if it is reahzed, compares with that due to the body 
forces on fluxoids which are discussed here. 



3 Effective Forces on Fluxoids 

3.1 Pinning Force 

The rotation of the neutron superfluid imphes that there exist angular-momentum-carrying 
neutron-vortex lines within the core of a neutron star. These are believed to form a regular 
lattice parallel to the rotation axis of the star which is in general tilted with respect to the mag- 
netic axis, ie. the direction of the fluxoids. Thus, a neutron vortex might pmsept a "pinning 
site" against a moving fluxoid, vice versa, should the two structures overlap llaCj'Ea; some 10"^^ 
fluxoids and 10^^ vortices should be, typically, present in a neutron star. In the steady-state, the 
n-vortices must be co-rotating with the charged component of the star, including the lattice of 
the p-fluxoids, at a rate Q. For a given rate f2 of spinning down of the star the radial velocity 
of the n-vortex outward motion, at the core boundary, would be 

Rc Rc^ , . 

where i?c(~ 9 x lO'"^ cm) is the radius of the stellar core, and il.s{^ ^) is the superfluid rotation 
rate. A superfluid normally spins down while maintaining a positive rotational lag uj{= ilg — > 
0) with its vortices. Hence an outward radial Magnus force .Fm = p^nRcLo would act on the 
vortices, per unit length at the core boundary, whejie ps is the neutron-superfluid density, and 
K — 2 X 10~'^cm^s~^ is the vorticity of a vortex lineEj. 

The Magnus force on n-vortices has been usually assumed to be balanced by the viscous 
forces being primarily caused by the scattering of electrons off the magnetized cores of these 
vortices in the interior of a neutron starc3. However in presence of the pinning forces it turns 
out that the latter would be the dominant force which will effectively equate with the Magnus 
force, resulting in a net zero force on the n-vortices. This force balance condition may be used 
to calculate the magnitude of the pinning force on n-vortices, as well as its direction which could 
be, in principle, in either (inward or outward) radial directions depending on the relative motion 
of the vortices and fluxoids. 

The effective pinning force per unit length of each fluxoid may be, in turn, evaluated by 
noting the equality of the mutual pinning forces between fluxoids and vortices at each pinning 
site, and then balancing the total pinning force communicated between the two lattices of vortices 
and fluxoids. However, at any given instant of time only a small fraction of the vortices would 
be directly interacting with the fluxoids. The remaining much greater fraction of them (of the 
order of the ratio of an inter-fluxoid spacing to the size of a pinning interaction region; see below) 
should reside in the inter-fluxoid spacings. Thus the total Magnus force acting continuously on 
all the vortices may, or may not, be communicated instantaneously to all the fluxoids, depending 
on the assumed rigidity of the n-vortices and other considerations. We omit a detailed discussion 
of the derivation of the pinning force on fluxoids for the different possibilities and write-down the 
final results. The pinning force Fn acting on a fluxoid, per unit length, is derived asE3, either 

F„ = "^Fu « 2<^o Ps Rc - 5.03^ dyncm-i (2) 

where — and = ^ are the number densities per unit cross section area of the vortices 

and the fluxoids, respectively, (pQ = 2 x 10~^G cm^ is the magnetic flux carried by a fluxoid. 
Be = lO^Bs is the strength of the core field in units of G, ui-q is the superfluid lag iv in units 
of lO^^rad s^^, and Ps is the spin period in units of s. Notice that the sign of lu determines the 
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sign of Fn for which, as well as for the other forces discussed below, the outward direction will 
be reckoned as the positive sense. Or, else, 

F^^^ ("^Fu) = 2.59 X 10-^^^ dyn cm'' (3) 

where d[ = 2.3 x lO^^Bg ^ cm is the inter-fluxoid spacing, and dp is the effective size of a 
pinning interaction region around each fluxoid. A value of dp = A* = 118 fm has been used for 
the assumed magnetic pinning mechanism (see below), where A* is the effective London length 
of the proton superconductor, being also a length scale for the spread of the magnetic field of a 
neutron vortex line. 



Critical lag: 

The magnitude of the force which could be exerted at each intersection by a vortex on a fluxoid, 
and vice versa, is limited by a maximum value /p corresponding to the given strength of the pin- 
ning energy Ep and the finite length scale of the interaction dp; namely Ep — /pdp. This implies 
a maximum limiting value that the Magnus force on the n-vortices could achieve, corresponding 
to a maximum critical lag Wcr which is given as 

Wcr = 1-59 X 10"^ Bg/^ rad s"^ (4) 

The critical lag is the magnitude of the lag when radial velocities of fluxoids and vortices are 
different; ie. uj = LUa- or lu — —lUci when the vortices move faster or slower than the fluxoids, 
respectively. However, during a co-moving phase when the force communicated between a vortex 
and a fluxoid at each pinning point is less than its maximum value, /p, the lag might have any 
value within the range —ujci < < t^ci - The pinning energy arising from the magnetic interaction 
of a fluxoid- vortex pair is estimated to be Ep ^^10"^ ergs. A different estimate for the pinning 
y, due to the proton density perturbation, gives a smaller value of i?p ~ 5 x 10~^ ergs 
However, both mechanisms result in similar values for the pinning force, hence similar lUct, 
since the interaction length, dp, for the latter (= ^p) is smaller than for the magnetic interaction 
(= Ap) by about the same ratio as the inverse of the pinning energies. 




3.2 Drag force 

An isolated fluxoid moving through the normal degenerate electron gas in the core of a neutron 
star is subject to the viscous drag force of the electrons scatteriiiE-off its magnetic field. The 
viscous drag force, per unit length of a fluxoid, is estimated to beua 

^" - - -7.30 X 10^ v- dyn cm^^ (5) 

64 Ep^Xp c 

where Vp is the velocity of the outward radial motion of the fluxoids in units of cm s^^, Uc ~ 
3. X lO^^cm"'^ is the number density of the electrons, and Ep^ = 88 MeV is the electron Fermi 
energy, corresponding to a total density p = 2 x 10^'^g cm~^, and a neutron number density 
71,1 ~ 1-7 X lO^^cm"'^ in the core. 

The expression for Fy in Eq. 5 is derived based on the assumption of independent motions 
for single fluxoids. However, for the typical conditions in the interior of a neutron star the lattice 
of fluxoids might be, as a whole (or least as bundles consisting of not less than ten million 
fluxoids), "frozen-in" the electron gasQ. This is because the mean distance between the successive 
magnetic scattering of electrons by fluxoids turns out to be many orders of magnitudes smaller 
than the mean free path of the electrons due to other events, and also that the deflection angle 
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at each scattering event is very small. A detailed treatment of the coherent electron scattering 
by the fluxoid lattice has been shown to indeed require an almost zero relative velocity between 
the electrons and the lattice. Therefore, the flux expulsion out of the core might be prohibited 
altogether except if electron-current loops across the core-crust boundary is realized. Uncertain- 
ties about the true distribution of the magnetic flux and the correct value of the conductivity in 
the crust, and also the possibility of a mechanical failure of the solid crust due to a build-up of 
the magnetic stresses, however, obscure any definite conclusion to be drawn Ej. 

Moreover, there are oth£i_reasons to suspect the above suggested frozen-in approximation 
for the fluxoids as a whole DO. For example, the finite volume of the fluxoid lattice and also 
the influence of the superconductor boundary effects on the motion of the fluxoids which have 
not been included in the above mentioned studies of the coherent scattering could as well have 
signiflcant new consequences. In addition, the motion of the incompressible electron fluid in the 
interior of a neutron star has been argued to be divergence freeQ. Any motion of the fluxoids along 
with the electrons, in the frozen-in approximation, must be, therefore, of the same (divergence- 
free) nature. This is however impossible for the uniform lattice of fluxoids during its outward 
motion since the lattice constant keeps changing. Hence, a compromise between the flux freezing 
and the divergence-free motion of the electrons has to be worked out if any flux expulsion is to 
be accounted for. 

Given the above uncertainties, as well as the lack of any other definite prescription for 
calculating the drag force due to the electron scattering, we use a value of as given by Eq. 5 
in our models. This choice is also supported by noticing that a tentative expulsion time scale 
derived for the case of coherent scattering, as implied by the Hall drift of the flux a t the base of 
the crust, turns out to be similar to that based on the single fluxoid approximationllalla. 



3.3 Buoyancy force 

The buoyancy force on fluxoids in a neutron star arises for reasons analogous to the case of 
macroscopic flux tubes in ordinary stars. Because flux tubes are in pressure equilibrium with 
their surrounding the excess magnetic pressure causes a deficit in the thermal pressure, and hence 
in the density, of the plasma inside a flux tube which make the tube to become buomot. The 
radially outward buoyancy force on a fluxoid, per unit length, can be expressed asHEJ 

/ (f>o yin(Ap/,fp) „ ^, , -1 
F^-[^ ^^=0.51 dyncm (6) 



3.4 Curvature force 

The tension of a vortex line (such as a fluxoid) implies that a curved geometry of the line would 
result in a restoring force, the curvature force, which tries to bring the line back to its minimum 
energy straight configuration. The concavely directed curvature force Fc, per unit length, on a 
vortex having a tension T and a curvature radius S is given as Fc = T/S. An outward moving 
fluxoid might be bent outward and becomes subject to an inward curvature force, since its end 
points are frozen in at the highly conductive bottom of the crust. Thus, for a fluxoid with a 
tension Tp = ((/)o/47rAp)^ ln(Ap/^p) the curvature force, per unit length, would be given as&a 

Rc f y ln(Ap/^p) _ i?c p A -1 

' ' -th — —0.35 dyn cm , (7) 



S \^Ti\J Rc S 

using a value of ~ ln2 for an assumed spatially uniform distribution of the fluxoids. 

The electron currents at the bottom of the crust are however subject to diffusion processes, 
the rate of which would set a maximum limiting speed Wmax for the fluxoids end points to move. 
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We, therefore, assume that whenever Vp < Wmax the fluxoids remain straight and no curvature 
force win be acting on them (Fc = 0), since their end points are also able to move with the same 
speed. In the opposite case, when Vp > Wmax, the fluxoids would be bent outward and the force 
Fc as in Eq. 7 would be effective. This maximum drift velocity t'max for the magnetic flux in the 
crust is estimated, based on the Ohmic diffusion alone, as 

W - - = 3.18 X 10-^ (t^] cm s-i (8) 



T 



V 10 V 



where R — 10^ cm is the radius of a neutron star, and t is the assumed time scale, in units of 
yr, for the decay of the magnetic field in the crust. A larger value for Umax may be expected if 
the Hall drift of the magnetic flux at the bottom of the crust is also taken into account. 

On the other hand, the collective rigiditv-of the fluxoids lattice due to mutual repulsive forces 
between them might react to a deformationll3. The force Fc associated with even a piece of the 
lattice of a size of an inter-vortex spacing (including some lO'' flux lines) would be, accordingly, 
so large that any bending of the lattice is effectively prohibited. The velocity of the fluxoids 
would then be constrained at all times by the condition Vp < Umax- We use alternative models 
to test this possibility as well. 

4 The Models 

The steady-state radial motion of a fluxoid, in the region of interest, is thus dcjiifpiiined from the 



balance equation for all the radial forces acting on it, per unit length, that isE^ 

Fn + F,+Fh + F, = (9) 

which may be rewritten in the form 

- ""'"^ Pvp,+S = (10) 



where parameters a, (3, and 5{= -Fb -I- -Fc) would have different values for the alternative models 
that we have considered, and Vp^ is the fluxoid velocity Vp in units of 10^^ cm s^^. Recall that 
w-e, which is the value of ui in units of 10~^ rad s~^, might have either positive or negative 
values, as is also the case with 5 in some of the models. 

This single equation includes two unknown variables uj and Up, and represents the azimuthal 
component of the Magnus equation of motion for the proton vortices. No radial Magnus force acts 
on the fluxoids, hence the right hand side is set equal to zero, because of the assumed co-rotation 
of the fluxoids with the proton superconductor. There exist however additional restrictions on 
the motion of the fluxoids which help to fix the value of one of the variables and solve Eq. 10 
for the other. Namely, for a co-moving state Vp = is given and lu could be determined. In 
contrast, when Wp(^ v^) is unknown uj is given as lo — ujcr or lu = —ujcr for Vp < or Vp > Vn, 
respectively. Furthermore, inspection of Eq. 10 indicates that it admits one and only one of the 
three different solutions, for the given values of v^, Be, and Ps at any time, namely 

LU = Lu{Vp — Vn) iff —l-^ci < UJ < LUcr 

Vp ^ Vp{uj ^ Wcr) iff Vp < Vn 

Vp = Up(cj = -LUcr) iff Vp > Vn 

The rate of the flux expulsion out of the core. Be ~ —^B^Vp, and the evolution of the stellar 

surface field Bg (with a decay rate Bg = — ) are hence uniquely determined from the above 

force balance equation, given the spin evolution of the star which determines (Eq. 1) the vortex 
velocity Vn at any time. 
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4-1 The results 



The computed time evolution of Vp and oj are shown in Fig. 1, together with Vn and Wcr, as is 
predicted in one of the models. Characteristically similar results as in Fig. 1 are obtained for the 
other tested models as well. The fluxoids motion in Fig. 1 is seen to follow three evolutionary 
phases in which they move slower, together, and faster than the vortices, successively. Transitions 




1 10 100 1000 10^ 10^ 10^ 10' 10^ 10^ 10'° 1 10 100 1000 10* 10^ 10^ 10' 10^ 10^ 10^° 

t ( ,r ) t ( yr ) 



Figure 1: The right panel shows the predicted time evolution of the lag co and its critical value uicv in a solitary 
neutron star according to one of our simulated models. The left panel shows the corresponding evolution of the 
velocities of the fluxoids Vp, and the vortices Vn. Initial values of Bg = 10^^'^ G, Be = 0.9Bb, and a value of 
T = 10^ yr have been used. 

between these successive evolutionary phases occur because of the reduction in (cx Clg) as well 
as the increase in Pg; a final co- moving phase of fluxoids and vortices might also occur for some 
choices of the initial conditions. 

The predicted evolution of the core and surface fields according to the same model is given 
in Fig. 2. A substantial decrease in the core field occurs at a time t> 10^ yr, which is expected 
for the typical average values of i'p< 10~^ cm s~^ during the earlier times, because ^ = ^ 
implies that a time period At ^ is needed for a major reduction in the core field to occur. 




1 10 100 1000 10* 10^ 10^ 10' 10^ 10^ 10'° 1 10 100 1000 10* 10^ 10* 10^ 10* 10* 10^° 

t ( yr ) t ( yr ) 

Figure 2: The predicted time evolution of the strength of the magnetic field in the core Be and at the surface 
Bb, and the spin period Ps in a solitary neutron star according to the same model as in Fig. 1. The right panel 
(a) is for an assumed value of r = 10^ yr and corresponds to the results in Fig. 1, while the left panel (b) is for 
T = 10^ yr. 

However, because of the very small magnitude of Vp (although > Vn) and also the reduced value 
of Be at later times Be does not change, substantially, afterwards. The surface field Bg responds 
to the change in Be on the assumed decay time scale r of the crust. The nontrivial role of the 
stellar crust in these field evolution models may be seen by comparing Fig. 2a with Fig. 2b, where 
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values of T = 10^, and lO'* yr have been used, respectively. A larger value of r tends to maintain 
the initial Bs, hence a larger Ps as well as a larger Vn, over a more extended period of time. 
Consequently, smaller final values of and Bs are predicted for the larger assumed values of r, 
as is seen in Fig. 2. 

Further consideration of the predicted time evolution of the forces on fluxoids, contrasted 
with the field evolution, reveals that the dominant "driving" force for the flux expulsion is the 
buoyancy force which is positive throughout the evolution. Accordingly, the overall role of the 
pinning force in the field decay of neutron stars turns out to be more like a "brake" , preventing 



the flux from being otherwise expelled too rapidly. Our conclusion about t. 
pinning force is, however, new and in contradiction with the earlier claims 



raking role of the 
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